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Abstract 

We examine the majorization properties of general thermal-like mixed states de- 
pending on a set of parameters. Sufficient conditions which ensure the increase 
in mixedness, and hence of any associated entropic form, when these parameters 
are varied, are identified. We then discuss those exhibiting a power law distribu- 
tion, showing that they can be characterized by two distinct mixing parameters, 
one associated with temperature and the other with the non-extensivity index q. 
Illustrative numerical results are also provided. 
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The rigorous concept of disorder derived from the theory of majorization 
P12][5]H] has recently received renewed attention in theoretical physics, par- 
ticularly in the field of quantum information [SIISIITIISII^ . The essential reason 
is that it is stronger than that based on standard entropic considerations. The 
basic idea is that a given probability distribution or density matrix can be 
said to be "more mixed" or "disordered" than another only when it is ma- 
jorized by the latter. This implies a higher entropy of the former, although the 
converse implication is not necessarily true. Majorization provides a natural 
partial ordering on probability distributions [2j and has consequently found 
many applications not only in mathematics, but also in other areas such as 
economy and computer science. Moreover, majorization relations are often 
naturally satisfied. For instance, in discrete classical systems, the joint distri- 
bution of two random variables is always majorized by the marginal distribu- 
tions, while in quantum mechanics, the global density matrix p of a separable 
(i.e., non-entangled) mixed state of a composite system is always more mixed 
than the local reduced densities pi of each subsystem [7], a property which 
can be violated by entangled states. These statements are stronger than the 
corresponding entropic inequalities [i.e., S{p) > S{pi)]. 
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The aim of this work is to examine along the previous hnes the majorization 
properties of general thermal-like mixed states depending on a set of parame- 
ters, discussing as application those characterized by a power-law distribution 
pU|llfl2j . The latter have in recent years been analyzed and successfully em- 
ployed in a wide range of contexts p^l4|15|16|17j . and can be derived within 
a generalized nonextensive thermodynamic formalism based on the Tsallis en- 
tropy [To]. One of the basic physical questions we want to answer is if such 
states do become more disordered, in the way determined by majorization, 
when the temperature (or some other fundamental parameter characterizing 
the distribution) is increased, as occurs with the standard Boltzmann-Gibbs 
thermal state. Such property would have far reaching consequences, in par- 
ticular that of ensuring a universal entropy increase, i.e., an increase in any 
consistent disorder measure, and not just in that employed in the construction 
of the state. It would also imply the increase of the expectation value of any 
increasing function of energy, and not just of the energy itself. Here we will 
prove that such property is indeed valid. 

For this purpose, we first identify the sufficient conditions that ensure the in- 
crease in mixedness of a general thermal-like mixed state when the parameters 
that characterize it are varied. We then show that states exhibiting a power- 
law distribution can be characterized by two distinct mixing parameters, one 
associated with temperature and the other with the non-additivity index q. We 
also discuss the majorization properties of escort distributions and the mixing 
conditions for generalized thermal-like states in the presence of constraints on 
non-commuting observables. Illustrative numerical results for a simple model 
are provided as well. Distributions with correct mixing properties can then 
be employed to investigate the effects of randomness in physical systems in a 
more complete way than that provided by the standard thermal state, allow- 
ing to explore different paths from the ground state to the full random state, 
all characterized by a proper disorder increase. 

We will consider a quantum system of finite dimension n. The eigenvalues pi 
of any density matrix p for such system {pi > 0, J2i=iPi = 1) ^ill be sorted 
in what follows in decreasing order {pi > pj if i < j). A density p is then said 
to be more mixed than a second density p' {p -< p') if the eigenvalues of p are 
majorized hy those of p' [2||3|H] : 

j j 

p ^ p' ^ Sj = J2Pi < s'j = Y.Pi ^ j = 1, . . . , n - 1 , (1) 

i=l i=l 

with Sn = s'j^ = 1. In such the probabilities pi are more "spread out" 

than the pj's, and can be written as a convex combination of permutations 
of the latter, i.e., p = J2aQaPa{p'), where Pa are permutations and > 0, 
J2a Qa = ^ . The state described by p is then more mixed or "random" 
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than that described by p' . Accordingly, the completely random state p = I /n 
(with I the identity) is more mixed than any density, while any density is 
more mixed than a pure state, i.e., I /n -< p -< |$)($| V normalized density p 
and pure state |$). It can be also shown that if p ~< p', p can be written as 
a convex combination of unitary transformations of p', i.e., p = J2aQaU}^p'Ua, 
with > and UlUa = I, and viceversa (Uhlmann's theorem [3]). If the 
dimensions of p and p' differ, the same definition ([1]) can be applied after 
completing with zeros the set of eigenvalues of the density of lowest dimension. 

Let us briefly discuss now the relation with entropy. Consider for instance the 
general entropic forms [T5|19|20|[^ 

n 

SJ{p)=TTf{p)=Y.f{P^), (2) 
i=l 

where / is a smooth strictly concave function {f'{pi) < f'ipj) if Vi > Pj) 
defined in the interval [0,1], satisfying /(O) = /(I) = 0. The von Neumann 
entropy S{p) = — Trplnp and the Tsallis generalization [TO] , 

S,{p) = TT{p-p'')/{q-l), g>0, (3) 

which approaches the von Neumann entropy for g ^ 1, are the most important 
examples. It can be shown that if p ^ p' ^ Sf{p) > Sf{p') for any f of the 
previous form [SP^] (the same holds for the Renyi entropy S*^ = ln[l + (1 — 
q)Sg{p)]/{l — q) [23], since it is an increasing function of Sg{p), as well as for 
any Schur concave function of p [4J). However, for a given /, the converse is not 
necessarily true, so that the concept of disorder implied by Eqs. ([1]) is stronger 
than that based on a particular choice of /. Nonetheless, the converse holds 
as follows: if Sf{p) > Sf{p') for any f of the previous form =^ p p' |3||22] . 
In other words, the hallmark of increasing mixedness is a universal entropy 
increase. Note, however, that Eqs. ([T]) define a partial order relationship, in 
the sense that given two densities p, p', it may happen that p -/< p' and p' 7^ p. 

Let us now consider a general mixed state p(A) depending on a continuous 
parameter A. We will say that A is a mixing 'parameter in a certain interval if 
p(A) becomes more mixed as A increases in this interval: 

p(A) < p(A') if A > A'. (4) 

This is equivalent, in the case of a smooth dependence, to the condition 
dsj/d\ < , j = 1, . . . , n — 1 , 

within this interval. The generalized entropy 5'/[p(A)] is then a non- decreasing 
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function of A for any concave /, as easily verified: 

'-^-'ffiilf'M-fin»)]>o, (5) 

since f'{pj) — /'(pj+i) < for / concave. Such states exhibit then an unam- 
biguous disorder increase for increasing A. 

As a general example, let us consider the escort densities [23j 

Pg = pyZ„ Z, = Ttp\ (6) 

associated with a density matrix p. It is easily seen that for g>0, A = l/gis 
a mixing parameter for pq, since the ensuing partial sums satisfy, for j < n, 

ds- j n 

l^f = E E HP^/PJ) < 0, A = l/g > 0, (7) 

^'^ i=i k=j+i 

where Pgi = pf/Zg are the eigenvalues of Pg {pgi > Pgj V g > if i < j). Any 
other decreasing function of q is of course a mixing parameter for pg as well. 
We have therefore -< Pg' if < g < q' for a given fixed density matrix p. 

Assume now that the system is described by a Hamiltonian H with energies 
z = 1, . . . , n, sorted in what follows in increasing order, and consider densities 
p(A) which satisfy the conditions: a) they commute with if, h) their eigenvalues 
are non-increasing functions of energy {pi > pj if Ei < Sj) and c) A is a mixing 
parameter for p(A) in a certain interval. In such another fundamental 

consequence of Eq. (jl]) is that not only the average energy {H) p = Trp(A)if, 
but also the expectation value of any non-decreasing function w of H {w{ei) < 
w{ej) if i < j), independent of A, is a nan- decreasing function of A: 



This automatically ensures a non-negative generalized "specific heat" c\ = 
d{H)p/d\ > 0. Let us also remark that if dsj/dX were positive for some j 
(and Pj > Pj+i, Sj < ^j+i), one could always find functions / and w of the 
previous forms such that Eqs. ([5]) and ([8]) become negative. In this way, one 
can in principle always witness the absence of proper mixing increase. 

We may now say that p(A) exhibits a thermal-like behavior if in addition, it 
approaches the ground state density /i/ni in some limit A — >• Aq (/i denotes 
the projector onto the ground state energy subspace and n\ its degeneracy) 
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and the state of maximum disorder //n in some other hmit A Aoo, with A a 
mixing parameter for Aq < A < Aqo- The most common example of a state of 
the previous form is, of course, the standard Bohzmann-Gibbs (BG) thermal 
state (we set in what follows Boltzmann constant k = 1) 

p{T)=exp[-H/T]/Z{T), Z{T) = Ti exp[-H/T] , T>0, (9) 



obtained from the minimization of {H)p — TS{p). It is well known that its von 
Neumann entropy S[p{T)] is an increasing function of temperature, which is 
usually taken as the basis for the statement that p(T) becomes more disordered 
as T increases. However, in the present framework this statement can be more 
rigorously formulated. It is easily shown that T is a proper mixing parameter 
for p(T) in the interval (0, oo), i.e., 

p(T) -< p{T') if T > r > 0, (10) 



as the sums of its first j eigenvalues satisfy, for j = 1, . . . , n — 1, 

^ = E E P^Pk{e^-eu)/T^<Q. (11) 
C'-' i=l k=j+l 



Hence, not only its von Neumann entropy, but also its generalized entropy 
S'/[p(T)], is an increasing function ofT for an?/ concave /. The average of any 
increasing function of energy is an increasing function of T as well. It can be 
shown that the generalized thermal density obtained from the minimization of 
{H)p—TSf{p), given by [20] Pf(T) = f'~^{H/T+a), where a is a normalization 
constant and the cutoff f'~^{u) = if m > /'(O) applies, also becomes more 
mixed as T increases [T7] . Eq. fllip remains valid replacing piPk by Zpip^ , with 
Pi = -[f"{pi)Z]-^ for Pi> and Z = - Ep,>o[/"(Pi)]~S being both positive 
for / concave. 

Let us now identify the conditions which ensure that A is a mixing parameter 
for a density of the more general form 

p(A) = g{H, \)/Z{\) , Z{\) = Tr giH, A) , (12) 



where g{e,X) is assumed to be an arbitrary smooth positive non-increasing 
function of e for e G depending on a parameter A, with [p{X),H] = 0. 

The variation rate of the associated partial sums can be shown to be, for j < n, 



9s- " ^ ^ ding 

-^ = E E PiPk[9x{£t) - hi^k)] , gx{£) = -Q^, (13) 
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with Pi = g{ei, \)/Z{\), which generahzes Eq. (fTTi) . A sufficient condition 
which ensures dsj/dX < for j = 1, ... ,n — 1 is, therefore, that g\{£) be a 
non- decreasing function of e for e G [£!,£„], i.e.. 

For instance, if 

p(A) = g{H/X)/Z{X) , (15) 

where A > and g is here a positive non-increasing smooth real function, 
g\{£) = —X^^u[\ng{u)]', with u = e/X, and Eq. (fT^ leads to the condition 

-{u[\ng{u)]y>0, (16) 

which is not necessarily valid (consider for instance g{u) = [1 + ln('u + 1)]^^ 
for M > 0). It is, of course, vahd in the BG case for A = T {g{u) = exp[—u]) as 
well as for g{u) = exp(— m'') Vr>OifM>0 {—{u[\ng{u)]'}' = r'^u^~^ > 0). 
Note also that if (?(0) > and g{oo) = 0, the state p(A) = g{H/ X)/Z{X), with 
H = H — Ell, will always approach the ground state density h/ni for A ^ 
and the fully mixed state I/n for A —>■ oo. It is, however, the more stringent 
condition f|T6l) which ensures that p(A) will in addition become monotonously 
more mixed as it evolves from the ground state to the fully mixed state. 

Note also that if A is a mixing parameter for the density ( fT2l) . it will remain 
a mixing parameter for the associated escort density ([6]) for g > 0, which 
corresponds to gg{H, A) = g{H, A)''. The sign of ^q^^ is left unchanged for A 
independent of q. Eq. ([7]) also follows from Eq. (fTBl) for g{H, A) g^Hy/^. 

As an important example of Eq. (|T2|) . we will examine the mixing properties 
of density operators characterized by a power-law distribution, which can be 
written in the form of the Tsallis distribution [TOfTT] 

p(g, T*) = [/ - (1 - q)H/T*]PlZ{q, T*) , (17) 

where T* > represents an effective temperature, q the non-extensivity index, 
H = H—Sil the energy measured from the ground state and = {u+\u\)/2. 
For g — > 1, p{q,T) approaches the BG distribution Eq. f[T7|) is obviously 
positive and fulfills previous conditions a) and 6) V T* > and g G 5R. Its 
eigenvalues Pi are strictly decreasing functions of energy for g > 1, but just 
non-increasing functions for g < 1 due to the cutoff that applies in this case 
{pi = if (1 - q)ei > T*, where = 5^ - 5i). 
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Fig. 1. Top left: The von Neumann entropy of the Tsallis distribution (jl7p (scaled 
to the maximum value S"^ = S{I/n)) for n = 100 equally spaced levels with spacing 
huj, showing that it is a non-decreasing function of both T* and q. Remaining panels: 
The scaled Tsallis entropy Sq/[p{q, T*)]/ in the same system, as a function of T* 
at fixed q (right panels) and as a function of q at fixed T* (bottom left panel), for 
different values of g', ranging from q' = 0.25 to q' = 2 in steps of 0.25. Dashed 
(dotted) lines correspond to g' > 1 {q' < 1), solid lines lo q' = 1 (von Neumann 
entropy). They are all increasing functions of T* and q. 

We will now show that both T* and q are proper mixing parameters for p(g, T*). 

1 

Defining g{e,q,T*) = [1 — (1 — q)e/T*]]r'' , we obtain 



doling 1 d'^lng e 

didf^ = (T* - (1 - q)ey - ° ' 'didi ^ (T* - (1 - q)ey " °' ^^^^ 

for £ > and (1 — q)s/T* < 1, so that according to Eq. f|T^ . p(g, T*) becomes 
more mixed as either T* or q increases: 



p{q,T*)^p{q,T'*)iiT*>T'*>0, p{q,T*) ^ p{q' ,T*) if q > q' . (19) 

The role of T* and q as proper mixing parameters constitutes then another 
fundamental property of the distribution (fT7|) . Any entropy Sf[p{q,T*)] (in 
particular Sg'[p{q,T*)] for any q' > 0) is a non-decreasing function of both T* 
and q in any system, as illustrated in fig. 1 for a truncated harmonic oscillator. 

Let us remark that Eq. (fT7|) exhibits a proper thermal-like behavior with respect 
to both q andT*, since in addition it approaches the ground state density Ii/rii 
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both for T* — i> at fixed q and for g ^ — oo at fixed T*, and tlie random state 
I /n both for T* — oo at fixed q as well as for g — oo at fixed T*. Actually, 
due to the cutoff for g < 1, p{q,T*) = Ii/ui already for T* < (1 — q)A at 
fixed g < 1, and g < 1 — T*/A at fixed T*, where A is the lowest non-zero 
excitation energy. Note also that p{q,T*) — I /n behaves linearly with H for 
large T* {g{e,q,T*) ^ 1—e/T* for T* ^ |1 — g|£r) but logarithmically for large 
g {g{e, g, T*) ^ 1 - ln[l + qe/T*]/q for g > ln[l + qe/T*] and g > 1). 

Eq. f|T71) can also be characterized by other mixing parameters. For example, 
rewriting Eq. (fT7|) as ps(7,/i) = [1 — sH / pY/'^ /Z{'y, fi), with 7 = |1 — g| > 0, 
/i = T*/|l — g| > and s = Sign(l — g), feot/i /i and 7 are as well independent 
mixing parameters for both signs of s (despite the increase of 7 with decreasing 



g for g < 1), since = [^{p - 36)^' > 0, = [l\p - se)]-' > 



for se/fi < 1. Thus, Ps{l,p) -< Ps{l,p') H p > p' and Ps{l,p) -< Ps{l',p) 
if 7 > 7'. These properties also follow from Eq. flTBl) in the case of p (for 

g{u) = [1 — suY-l'^) and from Eqs. (l6])-(I7j) in the case of 7. Again, we obtain 
a thermal-like behavior with respect to both 7 and p for s = ±1, with the 
ground state density approached for yU ^ or 7 ^ and the full random 
state approached for /i — * 00 or 7 — > 00 (provided yU > if s = 1). 

The actual thermal state derived from the non-extensive thermodynamic for- 
malism based on the Tsallis entropy and the minimization of the free en- 



ergy Fq = {H)p^ — TSq{p) p!2][T5] . where g > and pq is the escort den- 

sity ([6]), is also of the form (fT7|) but with T* related to the actual T by 
T = [T* - (1 - q){H)p^]/Zq (following ref. [15]). For g > 1, T is a direct 
increasing function of T*, as in this case OT/dT* = {rjq — l)/[(g — > 

(with r] = 2{0)p^{0-^)p^ - 1 > 1 and O = / - (1 - q)H/T*), so that T will 
also be a proper mixing parameter. It will remain so for < g < 1 provided 
the absolute minimum of the free energy Fq at each T is considered |24J , as in 
this case the entropy (and hence T*) cannot decrease with increasing T. 

So far all previous expressions are applicable in both the quantum and the 
classical discrete case. Let us finally briefly examine the majorization prop- 
erties of density matrices constructed from constraints on two or more non- 
commuting observables. We may for instance consider two observables Hq, Hi, 
with [Hq, Hi] 7^ 0, and a density of the form 



where g{u) is a positive non-increasing function and A,^ > for u = 0, 1, 
generahzing Eq. (ITSl) . In this case, Eq. (fT3|) should be replaced, for j < n, by 




p(Ao, Ai) = g{Ho/\o + Hj\r)/Z{\, \i) , 



(20) 



dsj 



3 n 



dhigjuj) 



dX 



H PiPk[9i-9k], 9i 



{i\Hy\i 



), (21) 



2=1 k=j+l 



8 



where Ui denotes the eigenvalues (sorted in increasing order) of J2u=o,i Hy/ Xu 
and \i) the corresponding eigenstates (in case of degeneracy we assume in ( I2T!) 
Hy diagonal within each eigenspace). Hence, we can ensure that will be a 
mixing parameter for (I2U]) if does not decrease for increasing values of Ui. 

For a single observable Hq, g^ = —X^^Uig' {ui) / g{ui) and the previous condition 
reduces to Eq. flTB]) . This is also the case when {i\Hy\i) is proportional to Ui 
{{i\Hy\i) = aui, with a > and independent of i), as occurs in simple systems 
such as a harmonic oscillator H = [p'^/m + kx'^]/2. In this case {i\p'^\i) = mEi, 
{i\x'^\i) = Ei/k, with Ei = hjj{i + 1/2) the oscillator energies {u = ^k/m), 
so that both the mass m and the inverse oscillator strength k^^ are also 
mixing parameters in p = g{H/T)/Z{T), provided T is a mixing parameter 
too (for Ao = m, Ai = k~^, we would have g^ = —^X~^[uig'{ui)/g{ui)], with 
Ui = Ei/T). In the general case, however, {i\Hy\i)/ui may depend on i in 
a non-trivial way, so that the mixing properties of p will require a careful 
analysis of the behavior of g^ . 

In conclusion, we have applied the theory of majorization to identify the rig- 
orous sufficient mixing conditions, as well as their main physical implications, 
for general mixed states of the form f|T2l) . which are summarized in Eqs. f|T^ . 
(1161) . ([5]) and ([8]). As application, we have examined those characterized by 
a power law distribution, and shown that they can be expressed in terms of 
two fundamental mixing parameters, which can be taken as T* and q in the 
representation (fT7|) . We have in particular identified the role of q in (fTTl) as 
a rigorous mixing parameter. The actual thermal state derived in the Tsal- 
lis non-extensive thermodynamic formalism was also shown to become more 
mixed for increasing T (with the above remarks applying for < g < 1), as 
occurs with the standard BG thermal state, ensuring in particular a universal 
entropy increase, i.e., dSf[p{q,T)]/dT > for any concave /. These results 
strengthen thus the robustness of the generalized non-extensive thermody- 
namic formalism. We have also discussed the majorization properties of escort 
distributions [Eq. ([6])] and derived sufficient conditions for mixedness increase 
in the presence of non-commuting observables [Eq. fl2T]) ]. 

Majorization theory enables then to derive very general inequalities with deep 
implications by simple means. Generalized thermal-like distributions with 
proper mixing parameters may also help to provide a more complete descrip- 
tion of the behavior of a correlated quantum system with increasing random- 
ness, reveahng aspects which could be hidden in standard BG statistics [17J. 
For instance, any system possessing a limit temperature Tc in standard statis- 
tics for some property present in its ground state and absent in the vicinity 
of the completely random state, will also possess a limit value of the mix- 
ing parameter in a generalized thermal-like distribution. Knowledge of such 
boundaries (like a critical curve T*{q) in f|T7|) ) may provide a new perspective 
for the classification of order-disorder transitions or crossovers. 
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